Free-streaming approximation in early dynamics of relativistic heavy-ion collisions* 
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We investigate an approximation to early dynamics in relativistic heavy-ion collisions, where after 
formation the partons are free streaming and around the proper time of 1 fm/c undergo a sudden 
equilibration described in terms of the Landau matching condition. We discuss physical and formal 
aspects of this approach. In particular, we show that initial azimuthally asymmetric transverse 
flow develops for non-central collisions as a consequence of the sudden equilibration. Moreover, the 
energy-momentum tensor from the free-streaming stage matches very smoothly to the form used 
in the transverse hydrodynamics, whereas matching to isotropic hydrodynamics requires a more 
pronounced change in the energy-momentum tensor. After the hydrodynamic phase statistical 
hadronization is carried out with the help of THERMINATOR. The physical results for the transverse- 
momentum spectra, the elliptic-flow, and the Hanbury-Brown-Twiss correlation radii, including the 
ratio -Rout /-Rsidc as well as the dependence of the radii on the azimuthal angle (azIfBT), are properly 
described within our approach. The agreement is equally good for a purely hydrodynamic evolution 
started at an early proper time of 0.25 fm/c, or for the free streaming started at that time, followed 
by the sudden equilibration at r ~ 1 fm/c and then by perfect hydrodynamics. Thus, the inclusion 
of free streaming allows us to delay the start of hydrodynamics to more realistic times of the order 
of 1 fm/c. 
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I. INTRODUCTION 



The heavy-ion data collected in the experiments at the 
Relativistic Heavy-Ion Collider (RHIC) suggest that the 
matter produced in relativistic heavy-ion collisions equili- 
brates very fast (presumably within a fraction of 1 fm/c) 
and its subsequent behavior is very well described by the 
dynamics of a perfect fluid |il,[3,[l,E[i,i,[3,[i- The most 
common explanation of these features is the assumption 
that the formed matter is a strongly coupled quark-gluon 
plasma (sQGP) 0. Another popular explanation as- 
sumes that the plasma is weakly interacting, however the 
plasma instabilities lead to a fast isotropization of mat- 
ter, which in turn helps to achieve equilibration [10]. In 
this scenario one argues that hydrodynamics is applica- 
ble already at the time when the system is isotropic, but 
not necessarily equilibrated [llj. Yet another explana- 
tion is based on the fully 3+1 dimensional parton cascade 
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that includes inelastic pQCD-based bremsstrahlung and 
its back reaction [l^, [13] ■ 

Differences in various theoretical approaches to the 
early-stage dynamics reflect our lack of precise knowl- 
edge concerning the mechanism of particle production 
and their early evolution. Moreover, our recent calcula- 
tion reproducing consistently and uniformly the hadronic 
spectra, the elliptic flow, and the HBT correlation radii 
[14| stresses the importance of the detailed shape of 
the initial condition from which hydrodynamics starts. 
Hence, the uncertainty concerning the dynamics seems 
to be intertwined with the uncertainty concerning the 
initial conditions. Clearly, both the early dynamics and 
the initial conditions should be eventually obtained from 
the early microscopic dynamics, such as e.g. the Color 
Glass Condensate theory [l^, [l^, [l3] ■ In practice, how- 
ever, the modeling of the partonic stage at the required 
precision is a very difficult task. 

The concept of early thermalization is especially in- 
triguing. In some cases, in order to obtain a consistent 
description of the particle spectra and femtoscopy, the 
hydrodynamic model requires the initialization time as 
short as 0.1 fm/c [l^. The commonly used argument for 
the early thermalization phenomenon is that the spatial 
eccentricity, resulting from the spatial transverse asym- 
metry of the colliding nuclei at non-zero impact param- 
eters and driving the formation of the elliptic flow, de- 
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creases with time, hence, the experimentaUy observed 
large values of V2 require an early onset of hydrodynam- 
ics and, consequently, an equilibrated state. In this paper 
we carefully reexamine this point of view. We investigate 
in detail an approximation to the early-stage dynamics 
in rclativistic heavy ion collisions consisting of the free- 
streaming (FS) of partons followed by a sudden equili- 
bration (SE) to a thermalized phase, which subsequently 
undergoes a hydrodynamic evolution. This FS-I-SE ap- 
proximation has been proposed by Kolb, Sollfrank, and 
Heinz fl^ . It has been further considered in an investiga- 
tion of the isotropization problem by Jas and Mrowczyn- 
ski _20J , as well as elaborated in the context of the early 
development of flow by Sinyukov, Gyulassy, Karpenko, 
and Nazarenko [III, [Hill. 

The FS-I-SE approach assumes that after the forma- 
tion stage the partons are first free streaming and later, 
around the proper time of 1 fm/c, undergo a sudden equi- 
libration described in terms of the Landau matching con- 
dition. We discuss the physical and formal aspects of this 
approach, which is the basic goal of this work. In particu- 
lar, we show that for non-central collisions, where the sys- 
tem develops spatial azimuthal anisotropy, an initial az- 
imuthally asymmetric transverse flow develops as a con- 
sequence of FS-I-SE. Moreover, we show that the energy- 
momentum tensor obtained from the free-streaming stage 
matches very smoothly to the form needed for the trans- 
verse hydrodynamics, where the longitudinal pressure 
vanishes [l^ . The inclusion of the partonic free stream- 
ing starting at the proper time of about 0.25 fm/c fol- 
lowed by the sudden equilibration allows us to delay the 
start of hydrodynamics to comfortable times of the or- 
der of 1 fm/c. In the calcultations presented in this pa- 
per we use the isotropic perfect hydrodynamics, as de- 
scribed in Ref. After the hydrodynamic phase the 
statistical hadronization 26] is carried out with the 
help of THERMINATOR The obtained physical results 
for the transverse-momentum spectra, the elliptic- flow, 
and the Hanbury-Brown-Twiss correlation radii, includ- 
ing the ratio i?out/-Rside as well as the dependence of 
the radii on the azimuthal angle (azHBT) [2^, are all 
properly described within our approach. Thus, the ap- 
proach consisting of FS-I-SE followed by hydrodynamics 
from r = 1 fm/c may be used to obtain the uniform de- 
scription of the soft hadronic data in a very similar way 
as in Ref. [14], where the hydrodynamic evolution starts 
right away at the early proper time of tq ~ 0.25 fm/c 
and no free streaming is present. 

The outline of the paper is as follows: In the next 
Section we introduce the basic concepts, the physical in- 
terpretation, and the kinematics of the FS-f SE approach. 
In Sect. Ill the structure of the energy-momentum tensor 
of free-streaming particles is analyzed in detail to show 
how the asymmetric flow is generated in this framework. 
The Landau matching condition is worked out in Sect. 
IV, while the physical results obtained with the hydrody- 
namic and statistical-hadronization codes are presented 
in Sect. V. We conclude in Sect. VI. Throughout the 
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partons equilibration hydro freeze-out hadrons 

FIG. 1: (Color online) Schematic history of a system formed 
at mid-rapidity in relativistic heavy-ion collisions, consist- 
ing of partonic free streaming, equilibration, hydrodynamics, 
freeze-out, and free streaming of hadrons to detectors. In the 
top panel the equilibration as well as freeze-out occur grad- 
ually, while in the approximate scheme depicted in the lower 
panel they occur suddenly. 

paper we use the units where c = h = 1 . 



II. FREE STREAMING FOLLOWED BY 
SUDDEN EQUILIBRATION 

A. Basic idea 

The basic idea to idealize the early stage of evolution of 
a system formed in heavy-ion collisions by a stage of coUi- 
sionless partonic free-streaming (FS) followed by sudden 
equilibration (SE), and then by hydrodynamics, has been 
proposed by Kolb, Sollfrank, and Heinz (l9| several years 
ago in the context of the development of azimuthally 
asymmetric flow. We visualize this approximation in Fig. 
[TJ The approach assumes a sudden but delayed transi- 
tion from a non-equilibrium initial state, consisting of 
free-streaming partons, to a fully thermalized fluid. Ever 
since it has been generally thought that the approach, 
which admittedly decreases the spatial asymmetry with 
time, leads automatically to a reduction of the elliptic 
flow, which hydrodymically develops from the azimuthal 
asymmetry of the density profile. However, the mecha- 
nism is subtle. While free streaming itself cannot gen- 
crate azimuthal asymmetry in the momentum distribu- 
tion, according to the common knowledge that interac- 
tions among produced particles are needed to achieve this 
goal, the sudden equilibration preceded by FS is in fact 
capable of developing azimuthally asymmetric flow. The 
point is that SE is a dynamical act, where the energy- 
momentum tensor of the system changes abruptly into 
a diagonal form (in the reference frame co-moving with 
the fluid element). That way space-flow velocity corre- 
lations are induced, which results in a collective elliptic 
flow, further enhanced by the subsequent hydrodynamic 
evolution. We discuss this crucial issue in a greater detail 
in the proceeding Sections, where we are equipped with 
the necessary formalism, in particular in Sect. ITV^ 

At this point we only wish to provide a qualitative 
argument for the development of flow in the FS-I-SE ap- 
proach. Consider an infinite slab shown in the left part 
of Fig. [21 which emits particles isotropically from sources 
denoted by dots. After some time particles reach the vol- 
ume element indicated by a box. For symmetry reasons 



3 




FIG. 2: (Color online) Schematic view of the development 
of azimuthally asymmetric flow from azimuthally asymmet- 
ric systems in the FS+SE approach. Arrows next to volume 
elements indicate the collective flow velocity. 



their average velocity, indicated by arrows, is perpendic- 
ular to the surface of the slab. Now, matching to per- 
fect hydrodynamics means that we are going to treat the 
fluid element collectively. In other words, one may imag- 
ine that the particles glue together and move collectively 
in the direction perpendicular to the slab. That way a 
correlation of position and flow velocity is generated. In 
a more realistic situation of the right-hand side of Fig. [2] 
one starts from an ellipsoidal distribution of sources. In 
that case one gets a larger flow along the direction of the 
shorter half-axis. A quantitative calculation is presented 
in the following sections. The point we wish to make 
here is the geometric origin of the azimuthally asymmet- 
ric flow in the FS-I-SE approximation, which reflects the 
original spatial asymmetry. As a result, for non-central 
collisions the starting condition of hydrodynamics, when 
delayed with FS-I-SE, includes the azimuthally asymmet- 
ric initial flow velocity. 

We remark that the inclusion of the non-zero trans- 
verse flow at the starting point of hydrodynamics is one 
of the possible ways to solve the RHIC HBT puzzle. This 
idea was first put forward in Ref. [2^ in the context of the 
thermal (hydro-inspired) models. Then it was discussed 
in Refs. [22, HI]. Quite recently, the importance of the 
initial flow has been strongly emphasized in Refs. (T8l.[30|. 



B. Physical interpretation 
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FIG. 3: (Color online) Straight-line trajectory of the parton 
from the proper time ro to the SE time r. 



decreases with time, or equivalently, the cross section in- 
creases. At first, it may seem quite paradoxical that a 
system which with time gets more dilute, becomes more 
and more likely to interact. However, recall that as the 
system gets more dilute, the average distance between 
partons grows, and as a result the strong coupling con- 
stant between colored objects increases. These confine- 
ment effects would make the partons in the system more 
likely to interact as the time goes on. Admittedly, as 
mentioned in the Introduction, the issue of thermaliza- 
tion is rather complicated, as it is difficult to assess if the 
initial gluon system, which is born very far from equilib- 
rium, has enough time to equilibrate before falling apart 
due to expansion. 



C. Kinematics 

Massl ess parto ns are formed at the initial proper time 
To = \Jt\ — and move along straight lines at the 
speed of light until the proper time when free stream- 
ing ends, r = \/t'^ — (cf. Fig[3]). We introduce the 
space-time rapidities ?7o = 5 log and r\ = \ log |^ . 

Elementary kinematics, following simply from the fact 
that the particles move along straight lines with the ve- 
locity of light, links the positions of a parton on the 
initial and final hypcrsurfaces and its four-momentum 
■p^ = (pTCOshy,pTCOS0,pTsin(/),p7'Sinhy), where Y and 
are the parton's rapidity and transverse momentum. 
We find 



The simplest interpretation of the FS+SE approach is 
to simply view it as an approximation to viscous hydro- 
dynamics. Indeed, instead of considering a complicated 
viscous system far from equilibrium, where microscopi- 
cally the scattering cross section of partons has a finite 
value, one employs an idealization, where initially the 
partons are free, and later develop an infinite cross sec- 
tion, which results in a sudden equilibration of the sys- 
tem. Certainly, the approach may work when viscosity 



Tsinh(77 — y) = rosinh(ryo — (1) 
X — xq ^ d cos 0, y — yo + d sin 0, 

d — = rcosh(F — rf) — \J + T^smh'^(Y — 77). 

The same equations are derived in plj through the use 
of the coUisionless Boltzmann equation. Due to Eqs. ([T]) 
the phase-space densities of partons at the proper times 
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As the ratio t/tq increases, the distribution ([6]) becomes 
more and more peaked, as shown in Fig. [D At large val- 
ues of t/to the rapidity distribution is so sharply peaked 
around Y = rj that effectively 



Fr^S{Y ~ri). 



(7) 



We note that this form is frequently assumed right away 
as an initial condition for subsequent evolution of the 
system. Here it effectively follows from the kinematics 
of free streaming and becomes better and better as r 
increases. 

With the form ^ Eq. Q yields 



FIG. 4: (Color online) The sharpening with time of the rapid- 
ity distribution of partons emitted at the space-time rapidity 
770, resulting from the free-streaming kinematics. The plot- 
ted functions are normalized to unity at the origin. At large 
proper times r the distribution approaches 5{Y — rj). 



To and T are related, 

d^N{T) 

dYcPpTdrjdxdy 



drjodxndyo 



d'^Niro) 



dY d'^PTdrjodxodyo 



S{rio —Y — arcsinh[ — sinh(77 — Y)]) x 
To 

S{x — xq — d cos (j))d{y ~ yo — d sin (p) . 



(2) 



d^N{i 



dYd'^pxdrjdxdy 



= n{x — At cos cj), 
5{Y - ii)f{pT). 



At sin < 



(8) 



where At = t — tq and /{pt) is a transverse momentum 
distribution of partons (note that d = At for Y — rj). 
In the calculations presented in the following sections 
t/to = 4. We read off from Fig. S] that with this ra- 
tio the spread in rapidity is a fraction of unity, hence 
very narrow and the approximation ([7]) is well justified. 



III. DEVELOPMENT OF ASYMMETRIC FLOW 



D. Narrowing of the rapidity distribution with 
time 

It is reasonable to assume for simplicity a factorized 
boost-invariant form of the initial distribution of partons, 



dY d'^pxdriodxodyo 



= n{xo,yo)F{Y ~r]o,pT), (3) 



where n is their density in the transverse plane. Below we 
will apply the profile obtained from the Glauber model 
as given by GLISSANDO ^31;], as well as a simple Gaussian 
profile 



n{xo,yo) = exp 



2a2 



262 



(4) 



where a and b depend on the centrality (impact param- 
eter) of the collision. 

When the emission profile F is focused near Y — rjo, 
for instance 



i^^exph(r-ryo)V(2a')], 



(5) 



with the rapidity width parameter a 1, then the kine- 
matic condition JT]) transforms it into a function of y — 77, 



F ~ exp ( — arcsinh^ 



— sin(y — rj) 
To 



/(2a2 



(6) 



A. Energy-momentum tensor from free streaming 

The energy-momentum tensor at the proper time t, 
rapidity 77, and transverse position (x, y) is given by the 
formula 



T^"" = / dYd^pT 



d'^Nir) 



dY d^PTdrjdxdy 



(9) 



= A 



d(j)n{x — At cos 4>,y — At sin ( 



/ cosh^?7 coshrjcoscj} coshrysin^ coshr^sinhTy \ 



coshyy cos ( 
cosh77 sin <; 



cos (p 
cos 6 sin ( 



cos cp sm ( 
sin^ (h 



cos (jysinhr] 
sin (/)sinh77 



y cosh?7sinh77 cos cfysmhrj sin ^sinhr; sinh rj ] 

where A is a constant coming from the pt integration, 
which factorizes out as a consequence of the approxima- 
tions adopted earlier. Due to the assumed boost invari- 
ance the further calculations may be carried out at 77 = 0, 
where we may drop the fourth row and column containing 
zeros, and write 

T'"' ^ A d(t)n{x ~ At cos (j),y - At sine/)) X 
Jo 




coscp 
cos^ 
cos 6 sin ( 



sm (p 
cos (f) sin ( 
sin^ d> 



(10) 
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FIG. 5: (Color online) Contour maps of the energy-density 
profile £ at the end of the free-streaming evolution, normalized 
to unity at the origin (top), and of the transverse velocity 
= a/i"! + (bottom). The initial profile is from Eq. Q 
with a = 1.9 fm and b = 2.5 fm, corresponding to centrality 
20-30% for the Au+Au collisions at RHIC. The initial and 
final proper times of free streaming are ro = 0.25 fm and 
r = 1 fm. 



B. Local rest frame 

Next, at each point {x, y) we pass to a local reference 
frame where the averaged three-momentum contained by 
particles in a volume element vanishes. The four- velocity 
of boost needed for the passage to this local rest frame 



is obtained from the condition 

T^^''{x, y)u,{x, y) = e{x., y)g^''u,ix, y), (11) 

with 

u'' ^il,v,,Vy,0)/Vl-v^, v^^vl+vl (12) 

Indeed, from the Lorentz covariance in the local rest 
frame Up^p(a;,?/) = (1,0,0) and Eq. (fTTj) takes the form 

T^lix.y) = e{x,y)g^^ (13) 

thus 

T™(x,2/)=e(a;,y), T^],(x,y) = 0, (i = 1, 2, 3). (14) 

Thus £ is simply the energy density of the system in the 
local rest frame. 

C. Energy-density and expansion velocity profiles 

We plot sample profiles of e and v for a non-central 
collision in Fig. [5l The plot corresponds to the proper 
time of r = 1 fm, with the starting proper time of free 
streaming at tq = 0.25 fm from the Gaussian profile 
of Eq. (HI). The width parameters are a = 1.9 fm and 
b = 2.5 fm, which corresponds to centrality 20-30% for 
the Au4-Au collisions at the highest RHIC energy. We 
note that both profiles are elongated along the y-axes. 
For the energy-density it reflects the shape of the initial 
density profile. For the velocity we also find a steeper rise 
along the x-axis than the y-axis, with clear anisotropy, 
or the space-velocity correlation, visible. 

The top panel of Fig. [5] shows the sections of the 
energy-density profile e, normalized to its value at the 
origin at time tq, i.e. Eq ~ £(0,0; r = tq), plotted as 
functions of the transverse radius p — x"^ + y"^ for sev- 
eral values of t. The solid (dashed) lines correspond to 
the in-plane (out-of-plane) directions. Certainly, the sys- 
tem is more elongated in the out-of-plane direction. As 
the proper time progresses, the system spreads out, and 
the value of e at the origin drops. The corresponding ve- 
locity of fiow is shown in the bottom part of Fig. [B] We 
note several features: the growth of the velocity with r, 
the nearly linear increase with p near the origin, and, im- 
portantly, a stronger flow in the in-plane direction. Thus 
the azimuthally asymmetric flow develops. This feature 
will be explained in Sect. IIIIP] below. 

We have repeated the above analysis for a Glauber- 
model profile as obtained from the GLISSANDO simula- 
tions [3l|. The applied model is the so-called mixed 
model, with 85.5% of wounded nucleons and 14.5% of 
binary collisions. The result for the Glauber profile is 
shown in Fig. [7l Comparing it to the Gaussian-profile 
case of Fig. [SI we note qualitative similarities, but also 
important differences. The flow near the origin develops 
more slowly in the Glauber result the drop of e 
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FIG. 6: (Color online) Sections of the energy-density profile 
£ (top) normalized to eo = e(0, 0;to) at the origin, and of 
the velocity profile v — + (bottom), cut along the 

X axis (solid lines) and ?/-axis (dashed lines) and plotted vs. 
p — x'^ + y^. The initial profile is from Eq. Q for centrality 
20-30% at TO = 0.25 fm. The e profiles are for r = ro = 0.25, 
1, 2, and 4 fm (from top to bottom), while the velocity profiles 
are for r = 1, 2, and 4 fm (from bottom to top). We note 
that the flow is azimuthally asymmetric and stronger along 
the X axis. 
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FIG. 7: (Color online) Same as Fig. [6] for the Glauber initial 
conditions from GLISSANDO. 



at the origin is also slower. We note that the two initial 
energy-density profiles we compare are normalized in the 
same way, as well as their second moments (the width 
parameters) are equal, i.e.. 



y dxodyon{xo,yo) = const. 
dxodyoxln{xo,yo) = a^, 
dxodyoyln{xo,yo) = 6^ 



(15) 



The difference shows at higher moments, with the 
Glauber profile being more flat at the origin than the 
Gaussian. As we have seen comparing Figs. [5] and [3 
these subtle differences make an impact on the flow pro- 
file and the strength of the flow velocity. These issues 
were emphasized in Ref. 14], where it was also shown 
that the obtained results, in particular i?out /.Rsidc , are 
better when the Gaussian profile is employed. 



D. Small time- and gradient expansion 

The qualitative features of the behavior presented 
above may be understood in terms of the low At and 
low p expansion. Expanding to lowest order in At, 

n(a; — Atcos(/), y — ATsint/)) = (16) 
n{x, y) — dxn{x, y) At cos — dyn{x, y) At sin 0, 

and integrating over 4> in Eq. (jlOp. yields the energy- 
momentum tensor in the form 

n 



-^^Ard^n \n | . (17) 

— \ATdyn \n 

The solution of Eq. pT|) gives to lowest order the eigen- 
vector u = (1, v), with the transverse velocity of the sim- 
ple form 

At Vn(x, y) 



v(x,2/) 



(18) 



3 n{x,y) 

For the Gaussian profile ([4]) this immediately results in 
the Hubble flow 

This behavior is clearly seen in Fig. [6] near the ori- 
gin. The validity of Eq. (|18p requires the condition 
At| Vn{x,y) \/n{x,y) ^ 1. 

IV. LANDAU MATCHING 

A. Matching to transverse and isotropic 
hydrodynamics 



As already stated in Sect. IIII B[ the boost of T'^" with 
the velocity found from Eq. pT|) yields T^p, i.e. the 
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energy density 



or to the perfect isotropic hydrodynamics, in which case 




X [fm] 

FIG. 8: (Color online) Same as Fig. [5]with the matrix Trf/s 
shown in a few points indicated by blobs. 



energy-momentum tensor in the local rest frame which 
satisfies the conditions (|14p . In Fig. [8] we have overlaid 
over the energy profile of Fig. [5] the explicit form of the 
matrix T^p/^ in a few points. We note that Tj^p is very 
close to the diagonal form 



^RF 




(20) 



In fact, for symmetry reasons Trf is diagonal along the 
X and y axes, and away from them it develops only small 
non-diagonal pieces. Also, the difference between T^^ 
and is small, at the level of a few percent. Interest- 
ingly, Eq. (PO)) has precisely the structure of the energy- 
momentum tensor of the perfect transverse hydrodynam- 
ics of massless particles [24!| , with the transverse pressure 
equal to e/2. Small departures from this form, present 
in our case, have the same structure as the shear tensor 
used to describe the viscosity effects in transverse hydro- 
dynamics [32] . We notice larger deviation in T^x and Tyy 
than in Txy, the same effect as in viscous hydrodynamics. 

The Landau matching condition amounts to replacing 
the free-streaming energy-momentum tensor to the form 
from perfect hydrodynamics. One may match to the per- 
fect transverse hydrodynamics, where we replace 



^RF 



/ 1 \ 

i 

i 

\ y 



(21) 



-'rf 



/ 1 \ 





\ 



(22) 



yO iy/ 



Admittedly, the matching (I2ip requires only a minor 
modification of the energy-momentum tensor, while (|22p 
employs a more significant change. This means that the 
free-streaming approximation may be rather smoothly 
linked to the perfect transverse hydrodynamics. 



B. Early generation of elliptic flow 

We interpret the Landau matching conditions ([2T|l or 
(|22p as a dynamical act. It is not a mathematical re- 
placement implementing an approximation. Rather, it 
is a simplified (recall Fig. [1]) description of the interac- 
tions among partons which occur instantaneously in the 
FS-I-SE approximation. The sudden equilibration causes 
the development of early elliptic flow. Following Ref. [l^ 
we consider the measure 



{Txx) + {Tyy) 



(23) 



The brackets indicate the space integration in the lab- 
oratory frame. Until SE occurs, Cp = 0, as without in- 
teractions the elliptic flow cannot develop, with the mo- 
mentum spectrum being an uncorrelated sum over the 
emitting sources. At the proper time Tmatch SE interac- 
tions occur, producing as a result the energy-momentum 
tensor of Eq. ([21]) or ([22]) . This act yields immediately 



a non-zero Sp. The subsequent hydrodynamic evolution 
may further increase the value of the elliptic flow coeffl- 
cient. The situation is depicted schematically in the top 
panel of Fig. [H 

In order to compute the value of generated by SE, 
we take T2 or T3 from Eq. dH]) or and for each 
volume element we go back from its local rest frame to 
the laboratory frame, which is necessary in order to apply 
the definition (|23|) . The result of this procedure is shown 
in the bottom panel of Fig.[9l We note that increasing the 
value of Tmatch rcsults in a larger flow coefficient. Thus 
FS-f SE does generate elliptic flow. At the same plot we 
also show the spatial eccentricity 



(24) 



which obviously decreases with time. From the view- 
point of hydrodynamics, this decrease of spatial asym- 
metry is compensated by the generated asymmetric flow 
from FS-I-SE. 
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FIG. 9: (Color online) Top: The schematic development of 
the partonic elliptic flow ep from the FS+SE approximation 
(top). Bottom: the value of the generated momentum asym- 
metry Ep plotted as a function of the matching proper time. 
The dotted (dashed) line corresponds to matching to isotropic 
(transverse) hydrodynamics, while the solid line shows the re- 
sults of hydrodynamics only, with no free streaming. The top 
curve shows the spatial asymmetry e. Same parameters as in 
Fig. El 
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FIG. 10: (Color online) The transverse-momentum spectra 
of pions, kaons and protons for c=20-30% (upper panel) and 
the elliptic flow coefRcient V2 for c=20-40% (lower panel). The 
darker (lighter) lines describe the model results for the case 
with (without) free streaming. Data from [ssl. l36l]. 



THE FOLLOW-UP EVOLUTION 



A. Hydrodynamics and THERMINATOR 

The energy-momentum tensor obtained with the 
FS-I-SE approximation is plugged into the hydrodynamic 
evolution. Here we use the perfect isotropic hydrodynam- 
ics, thus the matching condition (|22p is used. The equa- 
tion of state uses the lattice QCD simulations of Ref. [33| 
at high temperatures, T > 170 MeV, the hadronic gas at 
T < 170 MeV, and a smooth interpolation in the vicin- 
ity of 170 MeV. At the end of the hydrodynamic phase 
THERMINATOR simulations are carried out in order to im- 
plement the hadronic decays. Our scheme is described 
fully in Refs. 14, 34], so we do not provide any further 
details here. 



B. Physical results and comparison to data 

In the following we compare the results obtained 
with hydrodynamics only (starting at the proper time 
To = 0.25 fm) and the results obtained with free stream- 
ing from To to r = 1.0 fm, followed by SE and hydro- 
dynamics. In each case we start from the Gaussian pro- 
file Q. The results are shown in Figs. [TUl and [TT] with 



darker lines indicating the calculation with FS-I-SE, and 
the lighter lines with hydrodynamics only. We notice very 
similar results for the two considered cases, not to men- 
tion the very good description of the data. Larger free- 
streaming times (t — To 1.5 fm) spoil this agreement, 
as the flow becomes too strong. As described in Ref. [l3l , 
we have achieved a uniform agreement for soft physics at 
RHIC. In particular, the transverse-momentum spectra, 
the elliptic- flow, and the HBT correlation radii, including 
the notorious ratio i?out/^sido, are all properly described. 
The azimuthally-sensitive HBT correlations [stI . [sst are 
also correctly described within our framework. In Ref. 
[2^ we showed that our model calculations reproduce the 
full dependence of the HBT radii and their oscillations 
on the transverse-momentum and centrality. 

It is also interesting to look at the freeze-out hyper- 
surfaces for the two considered cases, shown in Fig. [T^l 
They have very much similar shape and size, which ex- 
plains again why the two schemes give essentially the 
same results. 

The practical observation following from our study of 
physical observables is that the inclusion of FS-I-SE may 
be used to delay the start of hydrodynamics. The phys- 
ical results are basically unaltered, since the dispersion 
of the density profile with time, resulting in milder hy- 



9 




0.2 0.25 0.3 0.35 0.4 0.45 0.5 
kj [GeV] 



FIG. 11; (Color online) The pion HBT radii R^idc , -Rout , 
i?iong, and the ratio Rout/Rside for central collisions. The 
darker (lighter) lines describe the results with (without) 
FS+SE. The data from H. 



drodynamic development of flow, is accompanied by the 
buildup of the initial asymmetric flow. 



VI. CONCLUSIONS 

Our analysis has focused on modeling of the early stage 
evolution within the FS+SE approach. We stress, how- 
ever, that our complete model (FS + SE + hydrodynam- 
ics + statistical hadronization) describes consistently the 
essential features of the soft hadron production at RHIC, 
including the pr-spectra, D2, and the pionic HBT radii, 
including their azimuthal asymmetry. The main reasons 
for obtaining such a good description, listed in Ref. p^ . 
were identified with the use of a realistic equation of state 
without the soft point, the Gaussian initial condition in- 
cluding the fluctuations of the initial eccentricity, as well 
as the inclusion of all known hadronic resonances in the 
statistical hadronization. This successful description of 
the RHIC data hints on a possible solution of the RHIC 
HBT puzzle 

In this paper we have analyzed the approximation for 
the early-stage dynamics of relativistic heavy-ion colli- 
sions which assumes the initial free-streaming of partons 
followed by the sudden equilibration (FS-I-SE approxima- 
tion). Our main findings are that such a modeling of the 
initial stage is compatible with the data describing soft 
hadron production at the highest RHIC energies. The 
decrease with time of the initial space asymmetry due to 
the free-streaming is compensated by the effects following 
from sudden equilibration, which lead to the formation 



of the radial and elliptic flow at the starting point for the 
hydrodynamic evolution. The net results of the FS-I-SE 
approach is that one may delay the initialization of the 
hydrodynamic evolution. Interestingly, unless the dura- 
tion of free streaming is not larger than about 1.5 fm/c, 
its specific value is irrelevant for the final physical ob- 
servables. Such insensitivity of the final results to the 
details of the initial evolution indicates that the assump- 
tion of the fast equilibration may be relaxed and replaced 
by a model where thermalization processes and building 
of the flow happen gradually. This issue of the persist- 
ing early thermalization puzzle might be resolved at a 
microscopic level only within an appropriate QCD-based 
kinetic model. 
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FIG. 12: (Color online) The freeze-out hypersurfaces for the 
calculation with hydrodynamics only (top) and with FS-I-SE 
followed by hydrodynamics (bottom). Same parameters as in 
Fig.El 
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